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Abstract 

We analyse the consequences of the usual assumption of a constant function to fit non-resonant 
decays from experimental Dalitz plot describing charmed meson decays. We first show, using the 

"> i^^Tr+Tr" decay channel as an example, how an inadequate extraction of the non-resonant 
contribution could yield incorrect measurements for the resonant channels. We analyse how the 
correct study of this decay will provide a test for the validity of factorization in D meson decays. 
Finally, we show how form factors could be extracted from non-resonant decays. We particularly 
discuss about the form factor that can be measured from the Df Tr^Tr+Tr"*" decay. We emphasize 
on its relevance for the study of the decay t — > v^Sn and the extraction of the ai meson width. 



1 Introduction 

Many body charm meson decays seem to be largely dominated by intermediate resonances. 
Experimental data have been studied using the powerful Dalitz plot technique which brings 
information on both the kinematics and the dynamics of the decay 

In a D meson three body decay, the intermediate resonant channels and the direct non- 
resonant one contribute to the final state. The Dalitz plot can thus present a complex interfer- 
ence of all these contributions. To extract them from the plot, one has to use appropriate fitting 
functions for each channel. Since the discovery of D mesons, data are fitted using a Breit- 
Wigner function for each resonance amplitude |0] while the non-resonant (NR) contribution 
has usually been considered as a phase space independent, constant function. 
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In a recent paper 0, we have shown that the last hypothesis cannot be safely considered for 
the study of D meson decays as it proceeds via a weak interaction. Indeed, in weak interactions 
at the partonic level helicity plays a central role; thus one could expect the amplitude of the 
reaction to have important variations within the phase space. In ref. 0, the NR contribution 
to — s> K^-K^Tx^ decay has been evaluated using factorization and an effective hamiltonian 
for the partonic interaction According to this calculation, the NR contribution does have 
significant variations along the phase space of the reaction. 

To extract data from the Dalitz plot, an adequate parametrization of the NR contribution 
is crucial. A correct extraction of the NR contribution could yield important information on 
the physics involving the decay; particularly, it can bring direct measurements of some form 
factors. Moreover, using an inadequate parametrization for the NR contribution, the whole 
decay pattern could be wrong: we could be ascribing to a given resonance those variations 
corresponding to the NR part. 

The purpose of this work is to present two examples concerning these ideas. First, we will 
analyse the decay K°7r+7r° , since its K*(892)°7r° partial decay width seems to be too 

large. Second, we will show how the ti'ti^ti^ decay is particularly well suited to extract 

a form factor which is relevant in r and ai-meson physics. 



2 The D+ K^7T+7T^ decay 

The resonant decay —>■ K* {892)'^7r~^ has been measured in two different ways, according 
to the detected final state: _ft'*°7r+) x B{K*° K°7r°), which is extracted from the 

Dalitz plot of the decay D+ K^n+n^; and B{D+ K*^7r+) x B{K*^ K-tt+), extracted 
from the decay —>■ K^tt^tt^. Marklll reported |^ an "apparent discrepancy" between the 
two measurements: B{D^ K*'^it^) = (5.9 ± 1.9 ±2.5)% when the final state is K^7r^7T~^ and 
B{D+ K*°7r+) = (1.8 ± 0.2 ± 1.0)% when the final state is K-tt+tt+. 

The last measurement has been confirmed by other experiments!^, |^, ^ while the decay 
— > A'^tt+tt" has only been measured by Marklll. It is then natural to think on a possible 
systematic error in the extraction of the A'* (892)*^ resonance from the — > A'^vr+Tr^ Dalitz 
plot. One possibility is that events coming from the NR contribution to the decay could have 
been incorrectly considered as originated from the ^*(892)'' resonant channel; thus, the latter 
has been artificially enhanced. 

To support this hypothesis, we are presenting here a calculation for the NR part of the 
decay —>■ K^n~^7r^ which shows up precisely an important bump near the K*{892)^ peak. 
The calculation is based in factorization |]T0| and in an effective Hamiltonian 0, |Tl| for the 



partonic interaction as in ref. 0. 

The effective partonic Hamiltonian is0, [11 



Heff = cos^ 9c[ai : {sc){ud) : +02 : {sd)(uc) :] (1) 
V2 
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where {qq') is a short-hand notation for g7^(l — 75)5'. The coefficients ai and 02 character- 
ize the contribution of the effective charged and neutral currents respectively, which include 
short- distance QCD effects. Figure 1 shows the six diagrams contributing to the amplitude 
^j^o^+^o- Using factorization one obtains the following decomposition [0] for the hadronic 
amplitude of the non-resonant decay 

-^m-.KO,o.+ = ^cosH,[a,{K'7T'\sc\D+){n+\ud\0) 
v2 

+a2(i^°|sc/|0)(7r°7r+|Mc|D+) + ai{K^\sc\D^) {7i\^\ud\0) 
+a2{K%^\sd\0){7T+\uc\D+)] . (2) 

Following ref. 0, the four contributions can be written as 

{K'Asc\D^){n^\ud\0) = -^/.m^Ff -°Ko.+ ,m^o.+ ) , (3) 

{7r^7r+\uc\D+){K^\sd\0) = ^fKml,Ff^\m\,^^,m\o^,) , (4) 
{K'\sc\D+){7f%+\ud\0) = -| { FhAmlo^^)fX4<o^+) (m|o.o - m^o.+ ) + 



FDXl"^7r07r+J/7r<'7r+l"^7r07r+J ^ ^ DK{^nOn+ ) 1^0^+ )\ " 2 i l^J 



and 



We have introduced the 3 invariants ^^11-0^+ = (p/fo -|- P-n+Y-, i^'ro^o = (P/?o + P-ko)^ and 
m^o^+ = (PttO +P7r+)^ and use has been made of the identity 

m|,+ + m^o + + m^o = m^o^+ + mj^o^o + mlo^+ . (7) 

The 1 / factors come from the wave functions. The ten form factors originate from the 
hadronic matrix elements|Q]; we will be back to them later on. 

Diagrams (a), (b) and (e) of Fig. 1 exhibit an external light pseudoscalar meson (P). This 
yields a contribution proportional to fprrip as one can see from eqs. @ and (|). The other 
diagrams, i.e., (c), (d) and (f), produce contributions proportional to m|) as one can see from 
eqs. (I) and (H), together with eq. (0). Thus the two ffist contributions in eq. (0) can be 
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safely neglected in favor of the last two. Moreover, the second term in eq. (|^) can be neglected 
as it is proportional to — ""^^o)- 

The eight form factors entering in eqs. (|^) and are written as 



(1 - q^lmlp j) 
where J = or 1 and P = K or -k, and 

/Aij(g')=/AB(0)(1 + AWM') (9) 

where z = + or and AB = Kit or tt+tt'^. According to the remark above, only six of them 
contribute to our calculation. Five of these form factors have been either measured from 
semileptonic decays, or calculated using lattice QCD or different quark models. There are no 
major discrepancies in the literature [^: F^^o(O) = 0.75 ± 0.1, m^^j^o = 2.0 ± 0.2; F^^(O) = 
0.75 ± 0.15, m\)^ = 2.1 ± 0.2; F^^{0) = 0.75 ± 0.15, = 2.2 ± 0.2; /|o^o(0) = 0.7 ± 0.1, 
Ato^o = 0.028 ± 0.002; /^o^o(O) = 0.7 ± 0.1, A^o^o = 0.004 ± 0.007. TheWth form factor 
entering in our calculation, f^+^o{q^), has neither been measured nor obtained using lattice 
calculations. One can only hint f^+^o{Q) calculating [|l^] the decay width 7r+ — n^e^Ve and 
comparing it with experiment, to get /^^o(0) = 1.4. 

Finally, the only measurement we have for the effective parameters ai and 02 come from 
the fit of two body charm meson decays. It has been found ai = 1.26 ± 0.10 and a2 = 
-0.51 ±0.10. 

We have performed a Monte Carlo simulation of the NR contribution to the decay — > 
K^n^TT^ using eqs. (H), (|]) and (|^). Figure 2 presents the Dalitz plot corresponding to the NR 
contribution to the decay, according to the calculation presented above. Figure 3 shows the 
density of events as a function of the invariant variable m|,o^o- It presents a pronounced bump 
centered in m'j^o^o ~ 0.65 GeV^. The figures shown have been obtained using the central values 
of all the parameters presented above. The bump remains unchanged doing any variation of 
these parameters within the region allowed by experiment and lattice calculation and a very 
large variations for the unknown parameters defining the f^+^o{q^) form factor. The bump is 
also unchanged within a large variation of the ratio 02/01. 

This strong stability of the bump in the simulation results shown in Fig. 3 is due to the 
following: it turns out that if the ratio 1 02/01! is not very large (in fact, smaller than about 
2.5), then the contribution of eq. (|^) largely drives the behavior of the ^|-o^o distribution. 
Thus, we can write 

M^+_Ko^o^+ oc F^;^(m^o^+)/+o^+(m^„^+) (m|o^o - m|o^+) (10) 

Since the form factors in ( PUD depend only on m^o^+, the "^l-o^o distribution of the events 
is thus almost independent of the various poorly known quantities associated with the decay. 
For comparison, we present in Figure 4 the m^Oj^o distribution of events when no dynamics is 
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assumed for the NR decay, i.e., Ai = const. The bump in Figure 3 is thus a robust signature 
of a factorization-based calculation. 

However, as some non-perturbative QCD effects - as final state interactions and soft gluon 
exchange - have been neglected, it is possible that factorization does not suffice to describe the 
decay. In the extreme case where the non-perturbative effects dominate the decay, the structure 
predicted above can be washed out. Thus, the experimental determination of a bump in the 
NR contribution to the decay ^°7r+7r° centered at mj^o^o ~ 0.65 GeV^ would be a test 

for the validity of factorization in D decays. 

The bump predicted by this calculation lies near the peak one expects for the Breit-Wigner 
distribution corresponding to the ^*(892)'^ resonance. If non-factorizable terms does not com- 
pletely eliminate the bump, many events originated from the NR decay K^7T~^7r^ have 
been probably incorrectly ascribed to the A'*(892)''7r"'" resonant channel. 

For completeness, we present in Figure 5 the Tnj^^^+ distribution of events for the NR part 
of the decay — > K^tt^tt^. We use the amplitude we obtained in ref. [Q. There is no bump 
near the the K*{892f squared mass. It looks more like a simple phase space distribution, as 
that of Figure 4. 

Thus, the difference between the NR contributions of the decays K^-k^-k^ and —>■ 

j^o^+^o could explain the different values reported for the decay K*{892)^tt^ according 

to the final state. It is an example of the indirect consequences of assuming inadequately a 
constant NR function to fit data. 



3 The Df 7T 7r+7r+ decay 

A correct extraction of the NR contribution to a given charmed many body decay could also 
have other important advantages. As we have shown above, the NR contribution to a given 
heavy meson decay is written in terms of various form factors. Thus, its correct extraction 
from the Dalitz plot could also be a way to measure those form factors within the whole phase 
space of the reaction. 

Two problems arise here. First, one has to accept that non-factorizable effects are small, so 
that the expression of the NR amplitude in terms of the form factors can be simply obtained 
using factorization hypothesis — as we did above and in ref. ||^. Second, even assuming 
the validity of factorization, those expressions are products of form factors and it is thus 
complicated to extract separately each of them. 

We present here an example in which these two problems are supposed to be not important. 
It is the case of the decay — > 7r^7r+7r+ . The main reason is that there is only one diagram 
contributing to the decay and it is an annihilation diagram. It is shown in Fig. 6. In this 



case, following Bjorken ideas[|T^], factorization looks natural. One thus expects the decay to 
be simply described by: 



= ^cos'eMOmDt){n^n^n-\A,\0) (11) 
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where 

m^iDt) = -ifoyo (12) 



and the second matrix element can be decomposed in four form factors|17]. The only remaining 



term is the axial spin one, i.e., similar to the one appearing in eqs. (|^) and One obtains 

G 



where m^_^, ^ (p.- +p .)^ and m^.^, ^ + p^,)\ 

The second problem raised above is naturally solved in this particular decay: the amplitude 
is proportional to just one form factor; thus, one can directly extract it from the plot. 

In two body decays of D mesons, amplitudes proportional to m|, only happen through 
spectator diagrams while those contributions coming from non spectator diagrams — as the 
one we are considering here — are proportional to the masses of the final state mesons, and 
thus less important. Since the amplitude of eq. (|13|) is proportional to m'j^^, in principle, it 
is not small. Nevertheless, if one assumes PCAC to be valid - due to the fact that final state 
quarks are light - one expects this decay to be small, and this can only happen if the form 
factor is negligible. However, the validity of PCAC in this context is not clear, as we will see 
in the following. 



The F4 form factor has never been measured and there are no clear theoretical predic- 
tions for it. Some authors]^ ^ proposed expressions based in models that are valid only 
for small values of the squared momentum transfered to the three pions, q^. However, in the 
decay 7r^7r+7r+, = m'jj^. 

The measurement of F4 will have important consequences on the understanding of r and ai 
meson physics. The ai width can be measured through the decay r z/^Stt but its value turns 
out to be 2 or even 3 times larger than the value extracted from other measurement s[^]. The 
value of the Oi width extracted from the decay r — > Ut-Sh strongly depends on the magnitude 
of a possible non-resonant decay, which is driven by F4, i.e. the same form factor involved in 
the NR decay 7r^7r+7r+. 

Experimental measurements!^ of the channel r — > UrSn cannot distinguish between mod- 



els predicting a large amount of PCAC breaking[22], i.e. a large F4, from those predicting 



a small amount of this breaking[g3[. Both a large and a small F4 are acceptable, but the 
extracted values of the ai width can vary by as much as a factor of two when fitting data using 
the first or the second kind of model. 

Thus, the correct extraction of the NR part of the decay 7r~7r"'"7r"'" can bring a first 

measurement of the form factor F4, clarifying the amount of PCAC breaking and then helping 
to extract the correct value of the Oi meson width. At present, the existent measurements of the 
decay 7r^7r^7r+ are not consistent: The branching ratio (BR) for the NR decay measured 

from the E691 experiment is 1.04 ±0.4 % - implying a large F4 - while preliminary results 
from E687p6l give a value of this BR as small as 0.121 ± 0.115 % - presenting a smaller F4. 



They have both been obtained using a constant function to fit the NR contribution. 
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4 Summary and conclusions 



In this work, we discuss on some of the consequences of our previous claim that NR con- 
tributions in D meson decays cannot be fitted with a constant as they usually are. We show 
here that important physics information is hidden in this contribution which has been loosely 
considered up to now. We present two examples to show the information one could obtain if 
the NR contribution were correctly extracted from the Dalitz plot. 

First, we argue that in the decay K^ti^tt^ events produced via the NR channel could 

have been assumed to be originated from the K*(892)° resonant contribution. Using a model 
based in factorization, we showed that NR have a bump near the K*{892)^ squared mass. This 
bump is very stable within a large variation of some poorly known quantities entering in the 
calculated amplitude. It is thus a strong prediction of factorization. 

Second, we claim here that an adequate extraction of the NR contribution from data could 
allow us to measure unknown form factors. The 7r^7r+7r+ is particularly interesting: its 

amplitude can be factorized to give a contribution proportional to the F4 form factor. This 
form factor drives the spin zero part of the axial current matrix element describing the decay 
in three pions. It is very relevant to the decay r z/^Svr. Different model predictions could 
be tested and the longstanding problem concerning the Oi meson width will benefit from this 
crucial information. 

Coming experiments on charmed meson decays are expected to measure — > K^tt^tt^ 
and D'^ —>■ 7r~7r+7r^ decays with high statistics. Using a non-constant function for the NR 
contribution when fitting the decay from its Dalitz plot, it will be possible to extract adequately 
the NR contributions. This will a) clarify the eventual discrepancy in the K*{892)^7t~^ 
decay width, b) test the validity of factorization technique when applied to D meson decays 
and c) bring a first measurement of the relevant form factor F4. 

One of us (RMG) would like to thank the Fermilab Theoretical Physics Department, where 
part of this work has been done. Two of us (IB and CG) want to thank the CNPq (Brazil) for 
financial support. 
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Figure 1: The six diagrams contributing to the decay D'^ — > K'^n^n'^ according to the effective 
Hamiltonian of equation (1). 



9 



^ 2.5 
t 2 



1.5 
1 

0.5 





I I I 



I I I 



I I I I 







0.5 



1.5 



m 



KtiO 



2.5 

GeV/c') 



3 

2\2 



Figure 2: The Dalitz plot for the NR decay D'^ X^tt+tt^. It has been obtained with Monte 
Carlo simulation weighted by \A4D+-^K0Tr+Tr°\^ ^ in equations (2), (5) and (6). 
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Figure 3: The mj^o^o density distribution for the NR decay — > ^^tt+tt". It has been 
obtained with Monte Carlo simulation weighted by \M.d+^ko-k+-ko\'^ ^ iii equations (2), (5) 
and (6). 
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Figure 4: Similar as Fig. 2, but for a flat decay, i.e., \M\'^ ~ const. 
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Figure 5: Similar as Fig. 2, but for the decay — > K 7r+7r+, using the model calculation 
developed in ref. [3]. 
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Figure 6: The annihilation diagram dominating the — > tt'tt'^tt'^ decay. 
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